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Abstract 

An upgraded (second) version of the package GENXICC (A Generator for Hadronic 
Production of the Double Heavy Baryons H cc , and Hj& by C.H. Chang, J.X. Wang 
and X.G. Wu, [its first version: in Comput. Phys. Commun. 177 (2007) 467-478]) is presented. 
Users, with this version being implemented in PYTHIA and a GNU C compiler, may simulate 
full events of the production in various experimental environments conveniently. In comparison 
with the previous version, in order to implement it in PYTHIA properly, a subprogram for the 
fragmentation of the produced double heavy diquark to the relevant baryon is complemented 
and the interphase of the generator to PYTHIA is changed accordingly. In the subprogram, 
with explanation, certain necessary assumptions (approximations) are made so as to conserve the 
momenta and the QCD 'color' flow for the fragmentation. 
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NEW VERSION PROGRAM SUMMARY 

Title of program : GENXICC2.0 

Program obtained from : CPC Program Library or the Institute of Theoretical Physics, 
Chinese Academy of Sciences, Beijing, P.R. China: www.itp.ac.cn/~zhangzx/ genxicc2.0. 

Reference to original program : GENXICC 

Reference in CPC : Comput. Phys. Commun. 177, 467(2007) 

Does the new version supersede the old program: No 

Computer : Any LINUX based on PC with FORTRAN 77 or FORTRAN 90 and GNU C 
compiler as well. 

Operating systems : LINUX. 

Programming language used : FORTRAN 77/90. 

Memory required to execute with typical data : About 2.0 MB. 

No. of bytes in distributed program, (including PYTHIA6.4) '■ About 1.5 MB. 

Distribution format : .tar.gz . 

Nature of physical problem : Hadronic production of double heavy baryons S cc , E bc and E bb . 

Method of solution : The code is based on NRQCD framework. With proper option, it can 
generate weighted and un-weighted events of hadronic double heavy baryon production. 
When the hadronizations of the produced jets and double heavy diquark are taken into 
account in the production, the upgraded version with proper interface to PYTHIA can well 
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generate the events in full. 



Restrictions on the complexity of the problem : The color flow, particularly, in the piece of 
programming the fragmentation from the produced colorful double heavy diquark into a 
relevant double heavy baryon, is treated carefully so as to implement it in PYTHIA properly. 

Reasons for new version : Responding to the feedback from users Q], we improve the gener- 
ator mainly by careful completing the 'final nonperturbative process' i.e. the formulation 
of the double heavy baryon from relevant intermediate diquark. In the present version, the 
information for fragmentation about momentum-flow and the-color flow, that is necessary 
for PYTHIA to generate full events, is retained although reasonable approximations are 
made. In comparison with the original version, the upgraded one can implement it in 
PYTHIA properly to do the full event simulation of the double heavy baryon production. 

Typical running time : It depends on which option is chosen to match PYTHIA when 
generating the events in full and also on which mechanism is chosen to generate the events. 
Typically, for the most complicated case with gluon-gluon fusion mechanism to generate 
the mixed events via the intermediate diquark in (cc)[ 3 Si]3 and (cc)[ 1 S'o]6 states, and to 
generate 1000 events, it takes about 20 hours on a 1.8 GHz Intel P4-processor machine 
if IDWTUP=1, whereas to generate 10 6 events it takes about 40 minutes only if IDWTUP=3. 

Keywords : Event generator; Hadronic production; double heavy baryons. 

Summary of the changes (improvements) : 1) We try to explain the treatment of the mo- 
mentum distribution of the process more clearly than the original version, and show how the 
final baryon is generated through the typical intermediate diquark precisely. 2) We present 
color flow of the involved processes precisely and the corresponding changes for the program 
are made. 3). The corresponding changes of the program are explained in the paper. 
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TABLE I: All considered mechanisms for step A, which are defined by the two parameters mgenxi 
and ixiccstate. Here the symbol gg-mechanism stands for the gluon-gluon fusion mechanism and 
etc.. 





mgenxi =1 


mgenxi =2 


mgenxi =3 


ixiccstate=l 


gg-mechanism, (cc)s( 3 Si) 


gg-mechanism, (60)3 ( 3 Si) 


gg-mechanism, (66)3 ( 3 <Si) 


ixiccstate=2 


gg-mechanism, (ccje^So) 


gg-mechanism, (6c)e( 1 S'o) 


gg-mechanism, (66)6 (^Sq) 


ixiccstate=3 


gc-mechanism, (cc)^( s Si) 


gg-mechanism, (bc)e( 3 Si) 




ixiccstate=4 


gc-mechanism, (cc^^So) 


gg-mechanism, (6c)g ( 1 Sq) 




ixiccstate=5 


cc-mechanism, (cc)g( 3 S'i) 






ixiccstate=6 


cc-mechanism, (cc^^So) 







I. MOMENTUM DISTRIBUTION OF THE PRODUCTION 

In fact, in the program we divide the production of a double heavy baryon H cc or S& c 
or into two steps: Step-A is up-to the production of a relevant double heavy diquark 
and Step-B is followed for the fragmentation of the double heavy diquark into the desired 
baryon. Therefore, in Step-A, there are three possible mechanisms: the gluon-gluon fusion 
mechanism (g + g), gluon-charm collision mechanism (g + c) and charm-charm collision 
mechanism (c+c), so in the program we need to fix one of them to produce the diquarks (cc), 
(be) and (bb) accordingly in term. All the three mechanisms and the calculation techniques 



for them are described in Refs.l2l, l3|, 14|. In Step-B, it is for the fragmentation of the double 
heavy diquark into the desired baryon. In the step we assume the intermediate diquark is to 
'decay' into the baryon plus soft parton(s) exclusively, e.g., either (QQ') pSijg — > ^qq' + Q or 
(QQ') [ 1( So]6 — > ^QQ' + <? + 9 (here Q, Q' denote c, 6-quark, q an light anti-quark, g a gluon). 
For convenience, in the program we name the mechanisms and intermediate diquarks in 
terms of the parameters as those in TABJH where mgenxi stands for the double-heavy 
diquarks, (cc) or (be) or (bb), and ixiccstate stands for the mechanisms. Note that in 
the previous version of GENXICC (we call it GENXICC1.0), we did not program how the 
diquark forms the relevant baryon, but alternatively we simply assumed that the relevant 
baryon is formed with 100% efficiency. 
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According to QCD confinement, the produced diquarks (QQ'), i.e. (cc), (6c), (bb), must 
be fragmented into relevant baryons by grabbing a light quark q (even suitable number 
of gluons g) with definite probability. Since the fragmentation for the heavy diquarks is 
absent from the available version of PYTHIA, thus in the upgraded version of GENXICC we 
program the fragmentation precisely and make its interphase still to suit PYTHIA properly. 
For consistency, in the upgraded version for the fragmentation we adopt the assumptions 
and method similar to those taken by PYTHIA [5] in the case for the fragmentation of a 
color-octet component (cc)§ into a colorless charmonium. Namely here the double heavy 
quark to grab a light quark (with gluons if necessary) from the 'environment' to form a 
colorless double heavy baryon with a relative possibility for various flavors of the light quark 
as«:G?:s:c~l:l: 0.3 : 1CT 11 . Hence in the program we introduce three new parameters 
ratiou (default=l), ratiod(default=l), ratios (default =0.3) so as to dictate the probability 
for a double heavy diquark to grab a light quark (antiquark) in forming the relevant baryon 
finally. These parameters may be changed by setting the values of the parameters in the 
parameter. F, when the relative possibilities for various flavors are assumed precisely. One 
more parameter nbound is naturally introduced to dictate which type of baryon: or 
Q+c (nbound=l), or Q bc (nbound=2), E bb ~ or Q bb (nbound=3) is to be generated 
from the relevant produced diquark. (nbound=4) is to derive the diquark results that can 
be generated by previous version (GENXICC1.0). The relative possibilities for the baryons 
S+, 5+ + and fi+ are decided by the value of ratiou, ratiod and ratios. More precisely, 
if the diquark (cc) [ 3 Si]g is produced, then it will fragment into S+ + with 43% probability, 
with 43% probability and £7+ with 14% probability accordingly, when default values of 
ratiou, ratiod and ratios are taken. 

Below we shall only take the hadronic production of S cc as an example to show how the 
generator GENXICC works, because the production for the baryon S;, c or H&& is similar. 

At the end of the Step-A, the final particles' momenta are set by using the phase_gen 
routine, that is based on RAMBO (Random Momentum Booster) program |f|, and the 
irrelevant phase space is integrated by VEGAS 7]. For the Step-B, we adopt the 'decay' 
method (as that to deal with the intermediate color octet (cc) states to produce J/ip in 
PYTHIA |5|]). According to the method, we start with assuming the diquark mass to be 
slightly bigger than that of the baryon (> 2m c ), so the diquark may 'decay' into a relevant 
baryon by emitting very soft partons (anti-quark and/or gluons), and the soft partons take 
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away very tiny momentum from the diquark. To keep gauge invariance of the hard process 
in Step-A, we must set (cc)-diquark mass to be 2m c exactly [3|, |4|, thus in the present case, 
we set the diquark mass to be 2m' c i.e. slightly bigger value, e.g. m' c > m c + m q /2 (m q is 
mass of the light quark), and let the mass of the produced baryon be less than but very 
close to 2m c + m q approximately. In the program, we introduce the parameters slqmass 
to stand for the soft-light-quark mass m g (q=u,d,s): qmassu for w-quark, qmassd for d- 
quark and qmasss for s-quark respectively. Since the Step-B in the case for the diquark 
(cc) pSijg being produced is of a (1 — > 2)-body process ((cc)[ 3 Si]3 — > E cc + q) and we have 
averaged the polarization of the produced diquark (cc)[ 3 .Sijg in Step-A, so we reasonably 
treat the fragmentation as that, for the double heavy baryon and the light anti-quark, the 
absolute values of the momenta is fixed by momentum conservation and the masses of the 
three relevant particles, whereas the direction of the momenta is isotropic in the diquark 
rest frame. It can be found that such a light anti-quark indeed affects very slightly in the 
momentum distribution of the generated double heavy baryon. We should note here that 
even though such a soft light quark gives slightly correction to the distribution in momentum, 
it must be involved so as to make the color flow of the process right for each event in full. 
If the color flow of the process were not correct, then the program would not be able to 
implement in PYTHIA to generate full events. As for the Step-B in the case for the diquark 
(cc)[ 1 5'o]6 being produced, that is of a (1 — > 3)-body process ((cc) [ 1 5'o]6 — ► 2 CC + q + g, one 
more gluon is needed) indeed due to the color-flow request as pointed out in the next section, 
we still treat it as a (1 — > 2)-body process for momentum distribution of the fragmentation 
by letting q and g always move together, "letting q and g always move together" is a very 
strong constraint but here q and g are very 'light' and carry the momenta very tiny, so it 
does not affect the distribution of the baryon, what we are interested in, substantially. 

To record all the useful data, we improve the data manipulation method. All the recorded 
data are still put in the directory data, while in this directory, we introduce 12 subdirectories; 
xiccu, xiccd, xiccs, xibcu, xibcd, xibcs, xibbu, xibbd, xibbs, ccdiq, bcdiq and bbdiq. These 
subdirectories record the needed data respective to the settings determined by the user 
in parameter. F, e.g. xiccu, xiccd and xiccs record data for f2+ s respectively; 

ccdiq, bcdiq and bbdiq record data for (cc)-diquark, (frc)-diquark, (frfe)-diquark respectively. 
We should note that the subdirectories ccdiq, bcdiq and bbdiq record the data just as those 
generated by GENXICC1.0. 
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II. THE COLOR FLOW OF THE PRODUCTION 



Ignoring the color flow of the fragmentation at all, still one can obtain incomplete infor- 
mation of the event for the baryon production, such as the observables: total cross section, 
transverse momentum p t and rapidity y differential distributions of the baryon etc. Thus, 
even GENXICC1.0, the previous version, can do all these jobs well. The present (upgraded) 
version deals with the color flow not only for the double heavy diquark production but also 
for the fragmentation carefully, so by complementing the upgraded version of the generator 
into PYTHIA properly, it can generate full events and record all information of the outgoing 
partons, besides that of the baryon itself. 

The color flow for S cc and H&& production is much more involved than that of "E^ produc- 
tion, so at the present, firstly we focus the color flow for S cc production precisely. Moreover 
to meet the needs in most cases, we only include the most important gluon-gluon fusion 
mechanism for the Sb c and production. The interesting readers can obtain the color 
flow for the S& c and S&& production by following the same method, but we note that for Sfe c 
production since the heavy quarks inside it are different flavored so its production via the 
diquark (be) with the quantum number as (be) ^So]:? or (&c)[ 3 Si]3 or (6c)[ 1 S , ]6 or (5c)[ 3 Si] 6 
(four possibilities) instead of (QQ)[ 3 Si]z or (QQ)[ 1 So\q (Q = c,b, two possibilities) for Eqq 
production. 



A. The color flow for (cc) [ 3 if?i]g being produced 

To show the color flow of S cc production, the mechanisms for (cc)[ 3 5i]g production and 
the followed fragmentation can be more precisely written as 

gg (cc) [ 3 5i]g + cc -»• (»+ + d) + cc or + u)+ce or (tt+ + sj + cc (1) 

gc -> (cc) [ 3 ^]3 + c -> (5+ + d) + c or + u) + c or (0+ + s) + c (2) 

cc -> (cc) fS^ + g^ (5+ +d)+g or + u) + g or + s) + g (3) 



According to the discussion in Ref. 8j, for the (cc)-diquark in 3 S\ state with the color 3, one 
of the two heavy c quarks emits a soft gluon and then this gluon splits into a light qq pair, 
and the (cc)-pair can combine the light quark q to form the baryon S cc (in \ccq) Fock state). 
To fit the needs of PYTHIA running, here we program the second step properly, i.e. to deal 
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FIG. 1: (Color on line) Typical color flow for gluon-gluon fusion mechanism (Left), gluon-charm col- 
lision mechanism (Middle) and charm-charm collision mechanism (Right), where the thick dashed 
line shows the corresponding intermediate diquark state (cc)[ 3 Si]3. Three colorful lines are for 
color flow lines according to PYTHIA naming rules. Three colorful lines are for color flow lines 
according to PYTHIA naming rules. The black lines are for Feynman diagrams. 

with the color flow of the fragmentation process carefully. 

As for the gluon-gluon fusion mechanism, there are six independent color factors for Step- 
A jjj: Cuj = (r a T b ^j G m jk, Cnj = 7^75 (T b T a ^j G m jk, C&j = ■^^{T a ) m j(T b ) n iG mn k, 

Giij — ^/fC^ )niG m nki G&j = T J^.Gmifc and Cftij = \^ ) m j^ m *^' 

where i, j = 1, 2, 3 are color indices of the two outgoing anti-quarks c and c respectively, and 
the indices a and b are color indices for gluon-1 and gluon-2 respectively. Here, the function 
G m jk either equals to the anti-symmetric e m jk when the (cc)-diquark is in 3 configuration 
or equals to the symmetric j m jk when the (cc)-diquark is in 6 configuration. While imple- 
menting Step-B, one may find that there are only two independent color flows for the case of 
(cc)[ 3 Si]3 with 50% probability each 1 . We draw the typical color flows for the gluon-gluon 
fusion mechanism, gluon-charm collision mechanism and charm-charm collision mechanism 
in FIG. ([I]) respectively, the other color flow can be obtained by gluon exchange. 
Now let us consider the three mechanisms in turn. 

For the gluon-gluon fusion mechanism as shown by FIG.(pJi), according to the naming 
rule of PYTHIA, there are three color flow lines: 

[0, 503] -> [502, 503] -> [501, 502] -> [501, 0] -> colorless bound state H+ c ++ /fi+, (4) 

1 This is slightly different from BCVEGPY, the generator for B c meson production, where there are five 
independent color factors in Step- A which lead to independent color flows accordingly [9(. 
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FIG. 2: (Color on line) Typical color flow for gluon-gluon fusion mechanism (Left), gluon-charm col- 
lision mechanism (Middle) and charm-charm collision mechanism (Right), where the thick dashed 
line shows the corresponding intermediate diquark state (cc)[ 1 S'o]6- Three colorful lines are for 
color flow lines according to PYTHIA naming rules. The black lines are for Feynman diagrams. 



[0, 504] -> [504, 0] -> colorless bound state H+ c ++ /fi+ , (5) 
[0, 505] -> [505, 0] -> colorless bound state H+ c ++/fi+ , (6) 

where the final sub-step in Step-B shows that the three partons c[501,0], c[505, 0] and 
g[504, 0] form the Fock state \ccq >. 

For the gluon-charm collision mechanism, the typical color flow as shown in FIG. ((lb ) is 

[502, 0] -> [501, 502] -> [501, 0] -> colorless bound state H+ c ++ /fi+, (7) 
[0, 504] -> [504, 0] -> colorless bound state S+| ++ /n+ , (8) 
[0, 505] -> [505, 0] -> colorless bound state S+| ++ /n+ . (9) 

For the charm-charm collision mechanism, the typical color flow as shown by FIG.flTb) is 

[501, 0] -> [502, 501] -»■ [502, 0] -> colorless bound state H+ c ++/fi+, (10) 
[503, 0] -> [503, 0] -> colorless bound state (11) 
[0, 504] -> [504, 0] -> colorless bound state S+' ++ /fl+ . (12) 



B. The color flow for (cc)f 1 5o]6 being produced 

In the case for (cc)[ 1 S , ]6 diquark being produced, except the color and the spin of the 
diquark, programming the Step-A is the same as that of (cc)[ 3 Si]3, but programming the 
Step-B is changed quite a lot: 
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• As pointed out in Ref . [8j] , the fragmentation of the color-sixtet diquark into the baryon 
is most likely via the Fock state \ccqg) instead via the Fock state \ccq) (the later need 
to change the heavy diquark spin, that it must be suppressed some). Note that as 
pointed out above, for the production of without 'identity particle rule', the role 
of the two Fock states \bcq) and \bcqg) played in the fragmentation is similar. So the 
difference in color flow for the \bcq) and \bcqg) does not cause different result in Eb c 
production. 

• Since the diquark (cc)[ 1 5 ]6 is in a color sixtet, in order to keep color flow of the 
production correctly, at least one more extra soft gluon, besides a light anti-quark, 
must be remained in the final state when the fragmentation is completed. It is because 
that, as the most simple gluon and an anti-quark can construct a color-sixtet: 

8(g)3 = 150603, (13) 

and the sixtet may fully 'absorb' the color of the diquark (cc) [ 1 5'o]6- Namely it is only 
due to the request to make color flow correctly that here an extra gluon is needed. 

• The most simple 'mechanism', which may be the biggest one, to generate the extra soft 
gluon is by emitting from another soft light quark or a soft gluon, which is necessary 
in the fragmentation sub-process (one typical example can be found in FIG.(J2])), so 
the extra so soft gluon causes negligible consequence in momentum distribution, even 
though all the subprocess are of non-perturbative nature. Hence in dealing with the 
momentum distribution among the baryon and jets etc, we can treat the soft anti- 
quark and soft gluon in the final state as a whole object 2 . Once more one may see 
that to introduce the soft anti-quarks and so soft gluons is only for keeping the color 
conserves (flows correctly) in the production. 

Some of typical color flows in the production via (cc)[ 1 S'o]6 are shown in FIG.(T3]), where 
FIG.(T5^i), FIG. ([2b) and FIG. ([2b) are for gluon-gluon fusion, gluon-charm collision and 
charm-charm collision mechanisms respectively. All possible colors for the production may 
be obtained by interchanging the connect points among the soft partons properly. 

2 It is reasonable, since numerically, one find that the soft quark itself gives negligible contributions to the 
momentum distributions for the rest particles in final state. 
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For the gluon-gluon fusion mechanism, i.e. the case as shown by FIG.fJJk), the three color 
flow lines are: 

a). [0, 503] -> [502, 503] -> [501, 502] -> [501, 0] 

— >■ colorless bound state Ej£; ++ /fi*, (14) 
6). [0, 504] -> [504, 505] -> [505, 506]-> [506,0] 4 

colorless bound state 3^ ++ /fi+, (15) 
c). [0, 507] [507, 0] ^ colorless bound state H+ c ++ /fi+, (16) 

where an extra soft gluon with color [505, 506] is introduced, in comparison with the case 
of (cc)[ 3 1S1J3, and connected to the second color flow line respectively. For the second color- 
flow line Eq.(15), both the gluon color g[505, 506] and the quark color g[506, 0] flow into the 
baryon. Thus here the equations Eqs. (14-16) correspond to that the four partons c[501,0], 
c[507, 0], g[506,0] and #[505,506] form the Fock state \ccqg > finally. It can be found that 
the color octet gluon g 3 [504, 505] and the color anti-triplet anti-quark g[0,507] can form a 
color 6 state so as to keep the color conservation of the process. 

For the gluon-charm collision mechanism, as shown by FIG. ([2b), we have 

a) . [502, 0] [501, 502] -> [501, 0] -> colorless bound state E+ c ++ /ft+ , (17) 

b) . [0, 504] [504, 505] -> [505, 506] [506,0] 4 

colorless bound state /ft^ c , (18) 

c) . [0, 507] -> [507, 0] -> colorless bound state H+ c ++ /n+. (19) 

The equations Eqs. (17-19) correspond to the four partons c[501,0], c[507, 0], g[506,0] and 
g [505, 506] form the Fock state \ccqg > finally. 

For the charm-charm collision mechanism, as shown by FIG. (J2b), we have 

a) . [501, 0] -» [502, 501] -» [502, 0] -> colorless bound state H+ c ++ /fi+ , (20) 

b) . [503, 0] -» [503, 504] -» [504, 505] ^ [505,0] 4 

colorless bound state E^T^/fi^., (21) 

c) . [0, 506] [506, 0] colorless bound state H+-++/Q+. (22) 

The equations Eqs. (20-22) correspond to the four partons c[502,0], c[506,0], g[505,0] and 
g[504, 505] form the Fock state \ccqg > finally. In Eqs. (15,18,21) the symbol " — > [• • •] — » " 
means a soft gluon is accompanied. 
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C. The color flow lines ended at a baryon and simplification with 'cheat' 

In the above two subsection, based on the color-flow decomposition method of Ref.[lo| 
we have designed a strict way to deal with the color flows for the baryon production. 

In fact since a baryon contains three valance quarks in different colors so as to form a 
colorless object, thus under the color-flow decomposition method, there must be three color 
flow lines being ended at a baryon. Namely it is different from the case of meson, where 
the color flow lines of the quark and anti-quark inside a meson are continued. To simplify 
the present situation, reminding the fact that 3® 3 = 603 in general QCD SU(3) color 
space: 1) Firstly, we can combine any two of the flow lines ended with two quarks into one 
anti-flow line ended with one anti-quark with a color 3 (the direction of the line is turned 
back completely) that is different from the two quarks (the third kind of color in respect 
those kinds of the two quarks); 2) Secondly, such anti-flow line obtained by the combination 
may be continued (connected) to the third color-flow line in the baryon; 3) Finally, as a 
consequence, the color-flow lines ended at a baryon become 'joined without ends' at all, 
which is the requirement of the color-singlet bound state. 

With the method for the color-flow lines ended at baryon, practically, the present version 
of PYTHIA (version 6.4) can generate full events with two or less independent color-flow 
lines in most cases well even a baryon being involved. 

When there are three or more independent color-flow lines in process, PYTHIA will stop 
at the step of running the 'final parton showers' or parton cascade radiations etc, and will 
present an error message to show that the color flow rearrangement is wrong during the 
parton's evolution process. For the present problem to generate the full events for the 
baryon production, it is found that in GENXICC all the mentioned production mechanisms: 
gluon-gluon fusion, gluon-charm collision and charm-charm collision, have exactly three 
independent color flows as shown by the last subsections when the color flow lines ended 
at the produced baryon are not treated as the above (to combine two of the three lines 
into one anti-color flow line). So with GENXICC1.0 being directly complemented, PYTHIA 
can not work well with these mechanisms at the step of 'final parton showers' or parton 
cascade radiations etc. To overcome the difficulty and to achieve full useful information for 
each event of the baryon production, the best approach perhaps is that PYTHIA itself is 
improved to deal with three or more color-flow lines' processes etc, e.g. to add the useful 
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intermediate diquark states with color 3 or color 6, such as (cc)[ 3 Si]3, (cc) pSo^ etc, into the 
particle tabular of PYTHIA, and properly to treat the color-flow lines ended at the baryon 
in fragmenting them into baryons accordingly. We suspect that it can be realized as done in 
a similar way as that in PYTHIA jlj] to treat color-octet mechanisms for the J/ty production 
without difficulty. 

Since GENXICC1.0, the previous version of the generator, without precisely strict color- 
flow lines being defined, cannot run well for the available versions of PYTHIA, and further- 
more, now there is no the improved PYTHIA as mentioned above, thus in GENXICC2.0 the 
upgraded version of the generator, we have to suggest another way so as to complete the full- 
event generation under a controlled approximate level (not to deviate from true too much) 
with the color flow lines ended at a baryon. Namely, we properly vary the color flows of the 
process as depicted in the last subsections, so as to make PYTHIA running well for final 
parton shower process. Aiming at the purpose, we precisely restrict the state constructed by 
either (ccq) or (ccqg) is equal to the desired colorless baryon, then I) to keep color conser- 
vation, all the outgoing particles' colors are arranged properly such that they form a color 
state as that of the incoming particles; II) The details in color of the constitute partons are 
ignored. PYTHIA is cheated to believe that the state constructed by either (ccq) or (ccqg) 
is colorless baryon, so we force it 'steady' without any evolution, thus any programs, which 
are called from PYTHIA, may run smoothly. In this way and only to introduce necessary 
soft enough light quark(s) anti-quark(s) and gluon(s), the momentum distribution and the 
other factors for the production will not be affected much. Moreover, we further simplify 
the program for the production in color as follows: 1) Ignoring the details about the color 
flow when the diquark fragments into a colorless baryon, the initial color-flow lines are re- 
duced to be less than or equal to two. More explicitly, for instance the three color-flow lines 
in Eqs.( Hll5ll6l) for the gluon-gluon fusion mechanism ((cc) pSijg) lead to two independent 
color-flow lines: [0,503] -> [502,503] -> [501,502] -> [501,0] -> colorless bound state 
and [0,504] — > [504,0] — > colorless bound state — > [0,504]; 2) A working assumption in 
program is made by setting the outgoing light anti-quark's color to be the same as that of 
the outgoing c 2 or (c) (see FIGfTJ in the case of gluon-gluon fusion or gluon-charm collision 
mechanism respectively; 3) In the production, to 'introduce' a soft anti-quarks and a soft 
gluon for the diquark (cc) ^Soje is only to keep the color conservation, so we may formally 
combine these two soft objects as an effective soft 'parton'. Therefore, as for the color 
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freedom-degree only, for the (cc)[ 1 S'o]6 we force the outgoing soft gluon and soft anti-quark 
to be combined into an effective one color-sixtet object (parton) formally, so the color flow 
for both (cc) [ 3 5i]g and (cc)[ 1 S'o]6 is treated in the same manner in the program. Namely 
in color, the effective soft 'parton' is q for (cc)[ 3 .Sijg and is (qg) for (cc) [ 1 5'o]6 accordingly. 
Moreover, in GENXICC2.0, we provide an interface of these soft objects to PYTHIA and 
users, who are interested in this part, can conveniently use these code to do the event evo- 
lution. We should note here that this simplification in treating the color flows cannot be 
applied to the charm-charm collision mechanism, i.e., we can not simultaneously simplify the 
three independent color-flow lines into two and also keep the process in color conservation. 
However, the contribution from this mechanism is small in comparison with the other two 
mechanisms at higher energy colliders as TEVATRON and LHC[3j fortunately, so in the 
generator we safely ignore the charm-charm collision mechanism at all. 

III. DISCUSSION 

To generate full events for the production of the double heavy baryons with GENXICC 
complemented into PYTHIA, as pointed out above, obviously the best way is to improve 
PITHIA with proper treatment the color-flow lines ended at the baryon. Whereas due 
to the fact that the improved version of PYTHIA is not available now, thus we 'have to 
suggest' another way with certain approximation, i.e. when writing down the program, the 
upgraded version has additionally made certain simplification with 'cheat' so as to generate 
the full events by applying GENXICC2.0, the upgraded generator, being complemented in 
PYTHIA, quite efficiently. However we should point out here that due to approximation 
and simplification with 'cheat' explained in the above subsection, the obtained information 
about the 'tiny jets', corresponding to the soft anti-quark and soft gluon(s) produced in 
fragmentation of double heavy diquark, may not be very reliable, although the information 
about hard jets and the baryon as well in the full events is reliable. 
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